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Q^' Abstract 

'^'h ' We use the density matrix formalisin in order to calculate the energy 

I/-S I level shifts, in second order on interaction, of an atom in the presence of a 

perfectly conducting wall in the dipole approximation. The thermal cor- 
rections are also examined when huJo/ksT = koXr 2> 1, where uq = koc 
^ , is the dominant transition frequency of the atom and At is the thermal 

CO ■ length. When the distance z between the atom and the wall is larger 

^ZJ ' than At we find the well known result obtained from Lifshitz's formula, 

whose leading term is proportional to temperature and is independent of 
c, h and feo- In the short distance limit, when 2 <C At, only very small 
corrections to the leading vacuum term occur. We also show, for all dis- 
/^^ . tance regimes, that the main thermal corrections are independent of fco 

■""•^ ' (dispersion is not important) and dependent of c, which means that there 

i-C I is not a non-retarded regime for the thermal contributions. 
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C ; 1 Introduction 

>Xh| In 1948, Casimir and Polder ^ considered for the first time the influence of 

retardation effects on the van der Waals forces between two atoms as well as on 
the force between an atom and a perfectly conducting wall. The non-retarded 
S^ . dispersive van der Walls force between two neutral, but polarizable atoms, had 

^ ' been explained previously by London in 1930 2'. Since Casimir and Polder's 

paper, forces between atoms or molecules and any kind of walls are usually 
called Casimir-Polder forces. These authors were motivated by experiments 
performed by Verwey and Overbeek [^ with colloidal suspensions which showed 
that London's interaction (that falls as 1/r^) was not correct for large distances, 
where the finite velocity of light should be taken into account. They showed 
that in the retarded limit the interaction falls as 1/r^ for two atoms (in contrast 
to London's result) and as l/r** for an atom and a perfectly conducting wall, in 
contrast with the short distance limit (that falls as 1/r^). Casimir and Polder 
obtained their results after lengthy calculations by using perturbative methods 
in QED ^. However, following a Bohr's suggestion, Casimir rederived the 
results obtained with Polder in a much simpler way, by computing the shift 
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in the electromagnetic zero-point energy caused by the presence of the atoms 
and the walls |5]- A month later, Casimir presented his seminal paper on the 
attraction between two parallel conducting plates which gave rise to the famous 
effect that since then bears his name [H]- 

In 1956, Lifshitz and collaborators developed a general theory of van der 
Waals forces [7). They derived a powerful expression for the force at finite 
temperature between two semi-infinite dispersive media characterized by well 
defined dielectric constants and separated by a slab of any other dispersive 
medium. They were able to derive and predict several results, like the variation 
of the thickness of thin superfluid helium films in a remarkable agreement with 
the experiments jH| ■ They also showed that the Casimir-Polder force is a limiting 
case when one of the media is sufficiently dilute such that the force between the 
slabs may be obtained by direct integration of a single atom- wall interaction |^ . 

Since then, many works have been done about this subject and the theory of 
van der Waals and Casimir-Polder interactions can be considered almost com- 
plete by now. Good reviews have been written on dispersive van der Waals 
interaction |1(J[ 1111 [T^ and many elaborated papers concerning level-shifts near 
surfaces have appeared, as for example, ^J^^EIEIj to mention just a few. 
It is worth mentioning that Casimir-Polder forces have been observed experi- 
mentally "^ . Higher multipole corrections JH| , roughness ^HI and corrugation 
of the surfaces [201 ' and the influence of the Casimir-Polder interaction on Bose- 
Einstein condensates 12] are some of the many branches of great activity on 
this subject nowadays, and the limiting cases of the Lifshitz formula are the 
starting point of the theoretical treatment of the majority of these problems. 

Here we propose an alternative way of computing the Casimir-Polder forces 
as well as van der Waals forces between single atoms which may be useful in the 
understanding of such interactions. The procedure to be presented is based on 
the density matrix formalism, from which we derive the master equation |221 for 
a small system interacting with a large one. It can be applied for any state of 
the field (not only vacuum or thermal states). Besides, it can take into account 
magnetic interactions (which are not considered in the Lifshitz formula) in a 
natural way. 

Consider a system 5-1-7^, where 5 is a small system and TZ a reservoir. 
Starting from an interacting Hamiltonian of the form V = — X)? ^j^j^ where 
Sj and Rj are S and TZ observables, the master equation comes from the pertur- 
bative expansion of the density matrix of the total system S -{-TZ. Going on up 
to second order, one can take the reduced trace over the Hilbert space of TZ and 
do the Markovian approximation. This approximation is due to the existence of 
two very different time scales, namely: the correlation time of the fluctuations 
of the TZ variables and the characteristic time evolution of 5, which is too large 
compared to the former. The final equation describes the evolution of the den- 
sity matrix of 5, made up of a free term plus a linear term that accounts for the 
coupling with TZ [2S| ■ One may then calculate the level shift of S and write it 
in terms of the linear susceptibility and symmetric correlation function of both 
small system and reservoir. This procedure permits us to distinguish the roles 
played by S and TZ to the interaction 22 . 

Here, the small system will be an electrically polarizable atom and the ra- 
diation field will be considered as the reservoir. In the dipole approximation, 
the Hamiltonian has the bilinear form required by the master equation. We will 
then calculate the interaction in both vacuum and thermal states of the field 



and investigate all possible distance regimes, defined by the main transition 
frequency of the atom and the thermal length (At = hc/kBT). 



2 The density matrix formalism 

In this section we will give a quick derivation of the master equation and obtain 
the general expressions for the level shifts. For details see references |22ll23l01) . 
Consider the time dependent perturbation: 



V{t)^-Y,S,{t)R,{t) 



(1) 



The equation for the evolution of the density matrix will be: 
d 



dt' 



p(t) = -Y^[S,{t)R,{t),p{t)] 



(2) 



Admitting that the density matrix oiS+TZ may be approximated by p = ps®PR, 
where ps is the reduced matrix of <S defined by ps = Tr^ [p] {pn is defined sim- 
ilarly) , integrating |(5Jl between to to t and making the adiabatic approximation 
(to —>■ — oo), one can show in first order on the interaction that the expectation 
value of an observable Sj at time t is given by: 



/oo 
dr x% (r) {Rk {t - r)) 
-OO 



Xl(r) = e(r)([5,(0),5.(-r)]} 



(3) 
(4) 



where Xjk (''") i^ ^^^ linear susceptibility of S respect to observables Sj and 
Sk EH and Q (r) is the Heaviside step function. The brackets mean that the 
ensemble average is over the non-perturbed states and (Sj) is the mean value 
of Sj at t = — cxo. In the deduction of last equations it was not necessary to 
approximate 7?. by a reservoir: the only assumption is that the total density 
matrix is a simple direct product of the density matrices of the subsystems. 

Considering now TZ as a, reservoir, one can choose {Rj) = for all relevant 
times without loss of generality, which causes the vanish of the first order term 
above. Another consequence for considering 7^ as a reservoir is that the second 
order correlations of TZ observables are very sharp in time (of order of the inverse 
of a cutoff frequency) compared to correlations of S observables (Markovian 
approximation). Then, we can write for the density matrix of 5: 
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yJabcd 



JabP^ab (t) + Y. ^<^bcdpil [t) 
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where <?,,(r) = Tr^^Jpfli?, (r) i?fe (0)], pf^ = {a\ps\b), Si, = {a\S,\b), 
hiOab ~ Ea — Ef) and \a) represents an eigenstate of energy of S with eigen- 
value Ea- Equation ^ is the master equation for S and gives all information 
about the system: the first term gives the free evolution and the second one 
accounts for the coupling with the reservoir. In order to obtain the level shifts 
it is necessary to consider only the non-diagonal terms which are coupled with 
themselves or, in other words, the terms in which the indexes are such that 
c = a and d = b. Then, from eq.lO, one may write: 

Jafcab = -Fab - zAab , /\^^-^p^^^V——————— (7) 

where Aqj, = Aq — Af,, Tab is a damping, V is the Cauchy's principal value, p^ 
is the statistical weight of the state |/i) which is an eigenstate of TZ with energy 
eigenvalue i?^ and V is defined by eq.lQJ^. 

Coming back to gjk (r), one may split it into its real and imaginary parts, 
which are related to the symmetric correlation functions and the linear suscep- 
tibilities. Then, we write A^ = A^J" + AJ^'' such that 



(8) 
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where SEa = hAa = hAl^ + TiA^ is the energy level shift of S in the state |a), 
x'^ (w) is the real part of the susceptibility of TZ and C^ (w) is its symmetric 
correlation function, C ^j" {ljo) and Xik"^ (^) ^re the equivalents of C^ {uo) and 
x'/ki^) for S in the state \a) and Rj^^ = {^\Rj\v). In the present formalism, 
equation (jS)) gives the contribution due to the polarization of the system under 
fluctuations of reservoir (/r) and equation (jSJ gives the reaction on the reservoir 
due to the fluctuations of the system (rr) . These equations will be the starting 
point for all calculations. 



3 Force between an atom and a wall 

Let us consider the standard Casimir setup: two perfectly conducting plates of 
side L separated by a distance ^, with L ^ i, and located at z = and z — £. 

^ Throughout this paper the greek indexes refer to reservoir energy eigenstates and the latin 
indexes refer to system energy eigenstates. 



Let us consider also that there exists an atom at a position z, with < z < ^. 
For a wave- vector k, the vector potential at the atom position is given by j26|: 



Ak||„(r||,z,i): 
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where uP' jc? = fc^ = fcj^ + [mi/i) , k., = k^x + fcj^y, r,. 



yij and n is a 



positive integer number. 

For an electrically polarizable atom the interaction Hamiltonian is: 

V{^,t) = -d(t)-E(x,i) 



(14) 



where d (i) is the dipole moment of the atom induced by the field. Since this 
Hamiltonian is bilinear in the atom and field operators, one may use expres- 
sions (jHJ and to calculate the interaction. Considering the atom as a two 
level isotropic system with transition frequency wq — ^oc, the correlations and 
susceptibilities of the atom and the field will be given by: 
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(15) 



kc + UJ fcc — cij 



2(nfc) + 1 



1- 



n^TT^ 2mT 



5{kc + uj) + 5{kc-uj) 



In the last equation, ao = 2|dp/3?ia;o is the static polarizability of the atom 
and (jifc) is the average number of photons for a given frequency and carries 
the information about the state of the field. Note that the susceptibility of the 
field is independent of its state and of h, which makes it a completely classical 
quantity (this happens under any boundary condition and in a linear medium) . 
Combining l|15(l and JHJ, summing over all possible modes and taking £ — > oo, 
one may write for the [rr) and (/r) contributions: 

5E" {z) = 5EI'' (z) = — / fc^a- (fc) G {2kz) dk ; 

T^ Jo 

he 



SEf"- (z) = SE^'' (z) + 5E{: (z) ; SE/^'' (z) = — / k^a+ (k) G {2kz) dk 

I" Jo 
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A first consequence of these results is that whenever the field is in the vacuum 
state, the (rr) contribution dominates at short distances and the (/r) contribu- 
tion dominates at large distances. For short distances {koz <^ 1), large values 
of k give the dominant contributions to the integrals: in this limit, a_ ~ 2/fc 
and a+ '^ — 2fco/fc^, which makes the (rr) term more important. At large dis- 
tances {koz ^ 1), only small values of k are significant: a_ ^ —2k/k'^ and 
a+ ^ 2/fco, and the roles of (/r) and (rr) contributions are interchanged. As 
another consequence, we should expect weak corrections to the vacuum term for 
short distances and strong changes in the behavior of the interaction at large 
distances if (uk) ^ 0, since the (/r) term is the only one that contributes. 

Focusing on the vacuum contribution, we may write for any regime: 

he kgao 
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(17) 
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Exact solution 
Non-retarded limit: 
k^z « 1 .0 

Retarded limit: 
kz» 1.0 
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Figure 1: Exact solution 1)17(1 and the asymptotic behaviors ((18|l and 119(1 for 
(koz) Vo{z), where Vo{z) is the vacuum contribution, in units of ftcaofcg. 



In the non-retarded limit {koz ^ 1), one may show that: 
Vo{z) = —-^ + 0{z ) 



(18) 



where too — koc. This result agrees with that obtained by the image method 0]. 
In the retarded limit (koz 3> 1), we have: 



^o(^) 



ao 



S-K Z^ ^ ' 



(19) 



which is the well known Casimir-Polder result PP . We show in Figure ^ the 
exact solution 1(17(1 and the asymptotic behaviors ((18(1 and 1(19(1 . 



3.1 Thermal corrections 



Let us now consider the radiation field in a thermal state at a temperature 
T . Recently, thermal effects on intermolecular forces at non-zero temperature 
have been considered in the literature. For instance, limits of validity of the 
London-van der Waals potential between two atoms (assuming dilute media) 
have benn established |27| . a general analytic expression valid at any tempera- 
ture T have been derived for the interaction of two polarizable electric dipoles 
PHj . or a Feinberg-Sucher theory [2^1 at finite temperature for two electrically- 
magnetically polarizable atoms in the non-dispersive regime have been devel- 
oped PO]. Here, we shall be concerned more specifically with the interaction 
between a polarizable atom and a perfectly conducting wall which is in thermal 
equilibrium at temperature T. In this situation, the atom interacts also with 
the corresponding thermal photons. 

As one can see from equations (|lfi|l . the contribution to the interaction is 
due only to the polarization of the atom by thermal fluctuations of the field: 

2nc f°° k^a+ (fc) 



Vt{z) 



= se!,: {z) - 



p(3hkc 1 



G{2kz)dk 



(20) 



In order to proceed it is convenient to introduce the quantity At — (3hc. It 
defines the only relevant length scale for the thermal behavior of the interaction 
when fco^T 3> 1 (this is a necessary requirement to make the approximation of 
the atom by a two level system a good one) in the same way that fco defines 
the only length scale for the vacuum contribution. Beyond At thermal effects 
become more and more important. Then, equation (|20|l may be written as: 
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where rjxo = koXr, xo = 2koZ and Ei (x) = J dt e ^ /t. For small values of 
zj \t, one can approximate the potential 1)21(1 by: 
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where N ^ rjxo is an integer and ( is the usual Riemann zeta function. The 
leading term of 1(22(1 is given by: 



Vr{z)^CiT)-^^^z' + 0(z*^ 



315 A^, 



(23) 



where C [T) ~ 1.38?icao/A^ is independent of z and, therefore, does not con- 
tribute to the interaction. This is a smaU repulsive ^ correction to the vacuum 
term (|Vt| ^ [z / \tY\Vq\) and gives a very good approximation for the true 
Vt[z) when z < 0.05At, as one can see in Figure O for koXr — 100, which 
is the usual value for optical transition frequencies at room temperature, this 
means that koz < 5; beyond this value, coincidently, the retarded effects become 
more relevant in the vacuum contribution (see Figure QJ. One may note also 
that (|23() is independent of ko and is c-dependent, which contrasts to (|18|l . This 
is an interesting result: at very short distances one could naively expect that 
the effects of the finiteness of the velocity of light might be completely ignored, 
which would imply an interaction independent of c. Though it is the case for the 
vacuum contribution, this does not occur for the thermal contribution: it is not 
possible to have an instantaneous field due to thermal fluctuations analogous to 
the Coulomb field at short distances. 
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Figure 2: Exact solution l|^ plus lfT7|l and the asymptotic Lifshitz behavior if^ 
for [z/\t) V{z,T), where V{z,T) is the total interaction, in units of ?icao/A^. 



For distances 2; ~ A^ or larger, expression l|22|) is not convenient. In this 
case, one can approximate the potential written in (|21() by: 
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^This result does not contradict the one found in reference I3UI for the thermal correction to 
the interaction between two atoms, Vr^-^ {r) = —CabT^/^+O (r) , since this last equation can 
not be used to derive 1231 by integration, because its validity is restricted to 1/fco <S r <C Ay, 
while 1231 takes into account contributions of all possible distances. 



where F'™' (a;) is the polygamma function of order m. Keeping the dominant 
terms of (|24|) and summing them with the vacuum contribution, one may write 
for the asymptotic total interaction: 



Viz,T) = Voiz) + VT{z) 



ksT ao 



(25) 



which is the weh known result obtained from Lifshitz formula [3 El ■ Looking at 
Figure El we see that the potential H25|l is in excellent agreement with the exact 
solution for z ~ At- For z > Xt the agreement is perfect (at room temperature 
this already occurs for z > 7.6fivr\). For comparison, one may note that the 
asymptotic behavior of the Casimir-Polder force given by H19|) fits the exact 
potential in a perfect way only for z > I.SAq, where Aq = 27r/fco. 
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Figure 3: Exact solution (|21|l . the short distance limiting case H23(l for < 
z/Xt < 0.05 and two like harmonic oscillator like potentials in the range 0.2 < 
z/Xt < 0.45 and 0.5 < z/Xt < 0.75 for Vt{z), which is the thermal correction 
to the interaction, in units of hcao/X'^. The vacuum term (|19|l is also shown. 



In the intermediate distance regime the thermal contribution to the interac- 
tion has a rich behavior. In Figure 13 we plot this contribution in the interval 
< z/Xt < 1-5. As one can see, there would be a stable equilibrium point 
(along the OZ direction) at z ~ 0.52At if this term were the only one con- 
tributing to the total interaction. For z < 0.52At, the thermal contribution 
to the interaction between the atom and the wall is repulsive, and in the other 
interval 0.52 < z/Xt < cxd, the thermal contribution has an attractive character. 
In the range 0.2 < z/Xt < 0.75 one may fit at least two like harmonic oscilator 
potentials: in 0.2 < z/Xt < 0.45 and 0.5 < z/Xt < 0.75. 

However, for the range 0.0 < z/Xt < 0.75, the vacuum contribution dom- 
inates completely over the thermal contribution: despite the complex thermal 
contribution behavior, this means that the total interaction is always attractive 



in this interval. For the other interval, 0.75 < z/At < cxd, the thermal term 
starts to compete with the vacuum term and becomes dominant from z ~ At to 
z = oo: but here, the thermal contribution is also attractive, so that the total 
potential T^(z,r) — Vq^z) + Vt{z) is again attractive. In conclusion, we can 
state that, for all distance regimes, the force between a polarizable atom and a 
perfectly conducting wall is always attractive, which is a well established result. 

4 Conclusions 

In this paper we discussed the standard example of Casimir-Polder force: an 
electrically polarizable atom interacting with a perfectly conducting wall. We 
used the expressions for level shift of the atom (considered as a two level system) 
due to the field (considered as a reservoir) calculated from the master equation. 
This approach allowed us to split the level shift into two terms: the reservoir 
reaction (rr) term which dominates at short distances and is independent of 
the state of the field; the fluctuation of the reservoir (/r) term which is more 
important at large distances and carries the information about the field state. 
We found the exact behavior of the interaction in the vacuum state in the dipole 
approximation and the correct limits at short and large distances (non-retarded 
and retarded limits respectively). 

We also considered the thermal state of the field. In this case we obtained 
the asymptotic behavior for distances z > At which gives a result in perfect 
agreement with that encountered in the literature |3 0] . We also found some 
results in the short distance limit {z < At) which are small corrections to the 
vacuum terms: \Vt\ ~ (z/At)''|Vo|. Considering only thermal corrections, we 
found a rich behavior for the interaction for z < At: repulsive and attractive 
zones with one stable equilibrium point at z c; 0.52At. In the short limit 
distance we found a thermal term depending on c, which is expected for fcoz ^ 1 
but not for koz <^ 1, where the velocity of light should be (naively) ignored: 
we showed that it is not possible to talk about an instantaneous interaction 
like Coulomb field for short distances when thermal corrections are taken into 
account, though this surprisingly occurs at sufficiently large distances (At ^ z). 
We think the formalism just presented may be useful for the computation of 
interatomic (intermolecular) interactions in many other situations, including 
magnetically polarizable atoms as well as other states of the radiation field. We 
hope to consider these problems in future works. 
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